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Abstract 

In the homogeneous and isotropic Friedmann-Robertson- Walker minisuper- 
space model, it is known that there are no Euclidean wormhole solutions in 
the pure gravity system. Here it is demonstrated explicitly that in Taub cos- 
mology, which is one of the simplest anisotropic cosmology models, wormhole 
solutions do exist in pure general relativity in both classical and quantum 
contexts. 
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Some time ago, the "Euclidean wormhole physics" stimulated an enormous excitement 
in the theoretical physics society due to the possible, provocative effects it may have on 
low-energy physics. For instance, the wormhole configuration may result in the generic loss 
of quantum coherence [2,3] or provide a mechanism for fixing the observable value of the 
cosmological constant to zero [1,2]. As for the quantitative aspect of this wormhole physics, 
the construction of explicit wormhole instanton solutions had been actively attempted [4- 
6]. To our disappointment, however, in the simplest minisuperspace model, homogeneous 
and isotropic wormhole solutions had been found only in the presence of somewhat exotic 
matter fields [8] but not in the pure gravity system. In this letter, we shall demonstrate 
that wormholes with the simplest anisotropic structure do exist in pure general relativity in 
both classical and quantum contexts. Classically, wormholes are Euclidean metrics which 
are solutions to the Euclidean classical field equations representing spacetimes consisting of 
two asymptotically Euclidean regions joined by a narrow tube. In the quantum sense, on 
the other hand, and particularly in the context of canonical quantum cosmology, "quantum 
wormholes" [8] may be identified with a state represented by a solution to the Wheeler- 
DeWitt (WD) equation satisfying a certain boundary condition [7] describing the typical 
wormhole configuration. Now consider the system of pure Einstein gravity with "negative" 
cosmological constant, i.e., the anti-de Sitter spacetime in 4-dim. Its Einstein- Hilbert action 
and the associated Einstein field equation are given respectively by 



where we introduced the Gibbons-Hawking surface term [9] on the boundary dM and A = 
— |A| < 0. We stress that the addition of negative cosmological constant here is not the 
necessary condition for the existence of wormhole solutions as we shall see in a moment. 
Then, in order to look for "anisotropic" (both classical and quantum) wormhole solutions, we 
choose to work in "Taub cosmology model" [10] (which is one of the simplest homogeneous, 
anisotropic cosmology models known) described by the metric 




(1) 
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ds'' = -N\t)de + e^-^'\e^^^%ja' (2) 

where the Euchdean signature metric can be obtained by replacing t by the Euchdean time 
r — it and thus —N'^{t)dt'^ by N'^{T)dr'^. N{t) and a{t) — e"^*) denote the lapse function and 
the scale factor respectively and Pij{t) = diag{(3, (3, —2/3) represents the degree of anisotropy. 
Thus this Taub model can be thought of as a special case of more general Bianchi type- 
IX model of anisotropic cosmology. Namely in the Bianchi type-IX model, the anisotropy 
measure Pij has two independent components /3+ and And the special case in which one 
sets —0 keeping = amounts to the reduction of the model to the Taub model. 
The non-holonomic basis {cr*} form a basis of 1-forms on a spacehke 3-sphere satisfying 
the SU(2) Maurer-Cartan structure equation 

da' = ^e^^ V^' A a'' (3) 

and can be represented in terms of 3-Euler angles 0<^<7r, O<0< 27r, < < 47r 
parametrizing as 

cr^ = cos ipd9 + sin ip sin 9d4>, 

— sin i/jd9 — cos t/j sin 9d(f), (4) 

— dip + cos 9d(l). 

The two independent metric functions, a{t) = e"^*) and P{t) thus can be regarded as two 
"minisuperspace variables" 7^ = (a, /3) and in terms of them the gravity action (dropping 
the Gibbons-Hawking surface term) can be written as 

dtNa'[^^{-{^r + P'} - {ly(/?) + ^A}] (5) 

where we defined the "curvature potential" V{P) as 

V{(3)^^-{e-'^-4e-'^) (6) 

and we redefined the action by multiplying it by an overall constant, S — > (37r/4G)S' and 
then rescaled the lapse function as N ^ 2N. For later use, we also provide, in the Lorentzian 
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signature, the expressions for the curvature scalar, it! of the 4-dim. space and the trace of 
the extrinsic curvature, K of a 3-dim. spacehke hypersurface ; 

Next, upon identifying the momenta tta = {Pa, Pp) conjugate to the minisuperspace variables 
7"^ = {a,P) as Pa — dS/dd — —ad/N and Pp — dS/d$ — a^P/N, one can go over to the 
Hamiltonian of the system via the Legendre transformation as 



S ^ j dtlADM = J dt{Pad + PpP - Hadm) 



with HADM^NHo + NiH\ Ni^O, (8) 
^0 = -H = la-'[-a'P'a + Pi + a'MP) + ^Aa'}] 

which (i.e., Hq) vanishes identically due to the time-reparametrization invariance in general 
relativity. This is called the classical Hamiltonian constraint, Hq = and, upon Dirac 
quantization, it turns into the Wheeler- De Witt equation to which we shall come back later 
on. In order to find classical wormhole solutions, which are Euclidean objects, from this 
point on, we need to work in Euclidean signature. The Euclidean gravity action can be 
obtained, via the analytic continuation, as 

where now the "prime" denotes the derivative with respect to the Euclidean time r while 
earher the "overdot" denoted that with respect to the Lorentzian time t. By varying this 
Euclidean action with respect to the lapse N{r), once again one gets the Hamiltonian con- 
straint 

^(^)=-^(f)= + [n« + iM^o (10) 

where (dp/dr) = —i{dp/dt) = —i{N/a^)Pfs using the Wick rotation rule, t = it while 
varying it with respect to minisuperspace variables a, one gets, respectively. 



1 .d^a. 2 .dp. 



1 ,d'^d. 3 ,da.,d3. 1 ,dV , „ 

Note, here, that the Hamiltonian constraint in eq.(lO) coincides with the rr-component of 
the Einstein equation, Rt-t- — ^gTT-R-\- Kg^-T — and to arrive at the expression for the eq.(ll), 
we made use of the Hamiltonian constraint in eq.(lO). Needless to say, constructing a general 
solution of the coupled equations (10), (11) and (12) is a highly non-trivial job. However, by 
restricting our interest to some particular class of wormhole solutions with specific character, 
we can have some insight into the existence of classical wormhole solutions. Namely, suppose 
we particularly look for wormhole solutions with very small/large anisotropy, 

V{P) = ^(e-^^ - Ae-^P) k = ^ 0), 0(/? ^ oo), (13) 

^ = li-e-^'^ + e-'P) - 0(/3 - or (5 ^ oo) (14) 

and with "constant" momentum associated with the anisotropy change, 

Pp = = const. (15) 

Imposing these conditions for the character of the solution, the eqs.(lO), (11) and (12) above 
reduce to, respectively, 

1 d ,N. 3 ,da, „ 

Next, we now fix the gauge associated with the time-reparametrization invariance as N{t) — 
1. Incidentally, then, the cq.(18) above is automatically satisfied and we are left just with 
the eq.(16) which is an ODE for the scale factor a(T) alone. (Eqs.(16) and (17), upon fixing 
the gauge N{t) — 1, become the same since one can get the latter by differentiating the 
former with respect to the Euclidean time r.) Thus what remains is to solve the eq.(16), 
i.e.. 
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for a(T) and then plug it into 



to obtain the behavior of anisotropy. Here, we set the "initial" time to be r — > — oo when 
yet neither the wormholes nor the baby universes were born and no anisotropy yet arises, 
i.e., oo) = 0. At this point, it seems worth noting that after fixing the gauge N{t) = 1, 
it becomes clear that the condition given in eq.(15), i.e., the constant momentum associated 
with the change in anisotropy, = const, corresponds to the condition for wormhole 
solutions whose anisotropy changing rate is inversely proportional to the scale (i.e., the size) 
of the wormhole, {dp/dr) — {—ijPp/a^, which seems quite a reasonable expectation. We 
now attempt to solve the eq.(19) for exact wormhole solutions, 
(i) wormhole solution with very small anisotropy {(3 — > 0) 

This case amounts to choosing k — —1, and then the eq.(19) admits an exact solution in 
the absence of the cosmological constant. Namely {da/dry — 1 — —Pp/a^ yields, upon 
integration, 

T = / dr'= / I = 21 

•^0 J<^) ^{a' + Pp){a'-Pp) 
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= y^F[arccos(^), - y2^^[arccos(^), + -^a^ - P^ 

where we chose that the wormhole throat, i.e., the minimum value of a(r) occurs for r = 0, 
namely (^)|^=o = and a{T — — ^JPp and F and E are the elliptic integrals of the 
1st and 2nd kind respectively. This solution can be identified with a classical wormhole 
configuration since it displays the asymptotic behavior a^(r) — > as r — >• ±oo. 
Now that we have found the exact wormhole solution in closed form. Then next, we discuss 
precisely what topology changing process this wormhole instanton solution describes. To 
this end, we note ; 
For T — > — oo, 
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Pi-oo) = 0, (f )|.^-oo = -^P,^^ - 

thus the curvature scalar there goes hke R{t — > — oo) — > 0, namely r —>■ — oo spacelike 
hypersurface is flat. 
For T = 0, 

a(0) = y^, (f)|.=o = 0, 
/?(0) = -zP,f_^ ^ = const., (f )|.=o = = 

thus the trace of the extrinsic curvature there is K{t = 0) = 3(^)|t-=o = 0, namely this 
T — spacelike hypersurface is a surface of vanishing extrinsic curvature and corresponds 
to the wormhole throat. And this surface has small, constant anisotropy and constant 
anisotropy changing rate. To conclude, this wormhole solution with very small anisotropy 
is an instanton configuration representing the topology change from the initial (r — > — oo) 
state with topology to the final (r = 0) state with the topology ® S^, i.e., flat space 
with an additional closed baby universe, 
(ii) wormhole solution with very large anisotropy (/3 — > oo) 

This case amounts to choosing k — and the eq.(19) admits an exact solution even in 
the presence of the cosmological constant. Namely {da/dry — 4|A|a^/3 = —P^/a'^, upon 
integration, 

^(^) , 1 , a'^da 



Jo Jw 

yields. 



(0) ^^Va^-r^' 



a(T) = [:^cosh(3«;T)]^/3 (22) 



K, 



where again we chose the wormhole neck to occur for r = and a(T = 0) = i'^iPpY^^ 
we defined = ||A| and r = Pp/n — PibJwm- 



We again discuss exactly what type of topology changing process this wormhole instanton 



solution describes. Thus we note ; 
For T — > — oo, 

a'ir) ^ ^ oo, (f ) 

Pi-Oo) = 0, (f )|._oo = -^^^^ - 

thus the curvature scalar there goes like R{t — > — oo) — >• 12k,'^ = 16|A|, namely r — > — oo 
spacelike hypersurface is the 3-dim. space of constant curvature with large radius, i.e., S^. 
For T = 0, 

P{Q) = -iPp ^ = const, {f)\r=o = -^P/3^ = 

thus the trace of the extrinsic curvature there is again K{t = 0) = 3(^)|t-=o = 0, namely this 
T — spacelike hypersurface is a surface of vanishing extrinsic curvature and hence is the 
wormhole throat. And this surface has large, constant anisotropy and constant anisotropy 
changing rate. Consequently, this wormhole solution with very large anisotropy is an in- 
stanton configuration describing the topology change from the initial (r — > — oo) state with 
topology to the final (r = 0) state with the topology ® S^, i.e., large 3-sphere with 
an additional closed baby universe. 

Having constructed wormhole instanton solutions, we now turn to their contribution to 
the topology-changing tunnelling amplitude. Namely, we evaluate the wormhole instanton 
action, I (instanton) , then the quantity, exp [—/(instanton)], would represent the semi- 
classical approximation to the topology-changing tunnelling amplitude. Thus by substitut- 
ing the (rr-component of) Einstein equation in (10) satisfied by the wormhole solution into 
the Euclidean action in eq.(l) or (9), one gets 

Hinstanton) = f'"' " (23) 



where we used {dp/dr) = -i{N/a?)Pi3, dr = da / N ^ \\A\a^ - V{(3)a'^ - P| and the 
Gibbons-Hawking surface term vanishes on each boundary, one at r = — oo and the other 



at T = as we have noted earlier. Note, first, that above expression for I{instanton) 
is independent of the lapse function N and hence the "physically observable" topology- 
changing tunnelling amplitude exp [—I{instanton)] possesses manifest gauge- invariance as- 
sociated with the time reparametrization invariance as it should. Unfortunately, it does not 
seem possible to carry out the integral in eq.(23) to obtain the precise value of the wormhole 
instanton action. Still, however, we can evaluate I{instanton) for wormhole solution with 
very small/large anisotropy we constructed earlier. Firstly, for the wormhole with very small 
anisotropy given in eq.(21). 



Hinstanton) ^ ^ ^^^^ H^f - f (25) 



where k and r are as defined earlier and Ai and A2 denote large cut-offs for asymptotic 
wormhole sizes a(— 00). In particular, the expression for the wormhole instanton action in 
eq.(25) indicates that the tunnelling amplitude for the topology change driven by the worm- 
hole solution in eq.(22) gets maximized for the smallest possible value of the cosmological 

constant, |A| = 3P|/4A^. 

We now turn to the discussion of anisotropic wormholes in the quantum regime. The formu- 
lation we shall employ for quantum treatment of wormholes [8] can be summarized as follows 
: we construct and study a minisuperspace model (again based on the Taub cosmology) of 
canonical quantum cosmology in which the main objective is to solve the Wheeler-DeWitt 
equation to find the universe wave function. Then to see if there is an excitation that can be 
interpreted as a "quantum wormhole" , we look for a particular solution to the WD equation 
with "wormhole boundary condition" that allows one to identify the universe wave function 
as representing an excitation corresponding to a wormhole state. And as a proposal for such 
wormhole boundary condition, we shall employ the one advocated by Hawking and Page 
[7]. According to them, wormhole boundary conditions can generally be classified into two 




(24) 



while secondly, for the wormhole with very large anisotropy given in eq.(22). 
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categories ; one for the "ground state" and the other for the "excited states" of quantum 
wormholes. Firstly, on the boundary condition for the ground state. When 0, the 

universe wave function should be regular to represent a non-singular 4-metric. And when 
\fh — > oo, the universe wave function should be damped, say, exponentially to represent an 
asymptotically Euclidean 4-metric, namely to represent the fact that there are no gravita- 
tional excitations asymptotically. Next, on the boundary condition for the excited states. 
When ^ ^ 0, again the universe wave function should be regular but it may oscillates for 
small- -\/^. And when \lh oo, the universe wave function should be damped. Thus now 
by applying the "Dirac quantization procedure" to general relativity, which is one of the 
most well-known constraint systems, the classical Hamiltonian constraint given in eq.(8), 
turns into its quantum version, namely the WD equation 

+ ^"1 - ^ + + 5A„^}l*[a,/,l = (26) 

where the suffix "p" represents the well-known ambiguity in "operator- ordering" and in 
passing from the classical Hamiltonian constraint to this quantum WD equation, we substi- 
tuted the conjugate momenta — —aa/N, Pp — a^(5/N by quantum momentum operators 
Pa — —id/da, Pp — —id/dp. And of course ^[a, /3] is the (physical) universe wave function. 
Again, we do not intend to challenge the general solution to this WD equation. Instead, we 
will be content with the universe wave function for quantum wormholes which have specific 
characters we considered earlier in the classical treatment. Namely, wormholes with very 
small/large anisotropy and with constant momentum associated with the anisotropy change 
described by conditions given in eqs.(13) (14) and (15). In this particular case, since V{(3) 
is replaced by a constant k (which is either —1 or 0), the WD equation above admits a sep- 
aration of variable. Namely by setting ^'[a, /3] = A{a)B{P), the WD equation decomposes 
into two sectors 

'^^+PlB = Q (28) 
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where we chose the constant involved in the separation of variable to be the conjugate 
momentum (which is a constant, too) since, as is manifest in eq.(28), it is the eigenvalue 
of the quantum momentum operator Pp = —id/d(5 which has to be the classical momentum 
Pp. Obviously, then, the solution of /3-sector is -B(/3) = [const. )e^^i^^ . Next, in order to have 
some insight into the behavior of the solution of a-sector, we view the a-sector of the WD 
equation given in eq.(27) as a Schrodinger-type equation with zero total energy Then the 
"potential" energy can be identified with 

U{a)^-a\k-UK\a^)-?l (29) 

which involves the term {—Pp/a^) representing an "abyss" in the small-a region. Since the 
total energy is zero, the emergence of this potential abyss reveals the fact that (part of) 
the universe wave function, A{a) should be a highly oscillating function of a in the small-a 
region. And this enormous oscillatory behavior for small scale factor a appears to signal 
the existence of "quantum wormholes" in the small-a region of the minisuperspace [8] and 
hence seems consistent with the existence of classical wormhole solutions we studied earlier. 
Thus in order to confirm this belief of ours, we now solve the a-sector of WD equation in 
eq.(27) exphcitly. 

(i) Quantum wormhole with very small anisotropy (/? — > 0) 

This case amounts to choosing k = —1 and then an exact solution to the a-sector of WD 
equation in (27) is available in the absence of the cosmological constant. 

(p + 1) d 9 Pa, . 

where Z^{x) denotes a Bessel function of order v — i^^AP^ — p^. For now Z,y{x) could 
be the Bessel function of the 1st kind or the 2nd kind (i.e., the Neumann function) or 
their linear combinations, i.e., the Hankel functions. However, imposing the appropriate 
"wormhole boundary condition" of Hawking and Page, it should be, for small 3-geometry 
(i.e., a — >• 0), the Bessel function of the 1st kind since 74(a) = a'^^J.i ^ 4p:2_p2 (|Q^) 
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a'^^a'^" — a~^P exp[i2\u\laa\ (where we used Ji,{x) — > j^yvX for x — > 0) with — 2P^ < 
p < for the regularity at a = 0, which indeed possesses highly oscillatory behavior for 
a — > 0. Next, for large 3-geometry (i.e., a — > oo), it should be the Hankel function since 
A{a) = a-5fiJ.(J^ ^__(ia2) ^ a-^^f+^^e-^"' (where we used li^J-\x) = J^{x)^iN^{x) 
x~2e*^ for X ^ oo), which indeed possesses rapidly damping behavior for a ^ oo as desired. 
Thus this solution *[a,/?] = A{a)e'^f'^ is a legitimate universe wave function of a quantum 
wormhole. 

(ii) Quantum wormhole with very large anisotropy (/3 — >^ oo) 

This case amounts to choosing A; = and an exact solution to eq.(27) is available even in 
the presence of the cosmological constant. 

Again, imposing the "wormhole boundary condition" of Hawking and Page, this general so- 
lution of the Bessel equation should be, for small 3-geometry (a ^ 0), the Bessel function of 
the 1st kind since A(a) — a~^^Xi^ ^p2_p2 (^|\/^Q^) a'^^a^" — a~5Pexp [i3|i/| Ina] with 
V — i^^APg — now and — 2P^ < p < for the regularity at a = 0, which exhibits highly 
oscillatory behavior for a — >■ 0. Next, for large 3-geometry (i.e., a oo), it should be the 



I 2 / |A| 3 

Hankel function since Aia) — a'^^H^]^ , — , — o(^iv^a^) a~^(p+3)e 3 " which ex- 

^ ' n\J^^l-^ 3 V 3 ^ 

hibits rapidly damping behavior for a — > oo as desired. And this solution ^^[0, /3] = A{a)e^^<^^ 
can be identified with a legitimate quantum wormhole wave function as well. Finally, we 
note that the wormhole wave functions given above in (i) and (ii) can be identified with 
the ones representing "excited states" of quantum wormholes we stated earlier as they are 
regular at a = and possess oscillatory behaviors for small-a. 

To summarize, in this letter, we demonstrated explicitly that Euclidean, anisotropic worm- 
hole solutions do exist even in pure Einstein gravity system. This, among others, strongly 
supports our intuitive expectation that there should be generic quantum fluctuations in 
spacetime alone even in the complete absence of matter. Noticing that homogeneous, 
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isotropic wormhole solutions do not exist in pure general relativity, one may naturally sus- 
pect that the anisotropy degree of freedom could play a decisive role in switching on the 
wormhole configuration. This is indeed the case as one can see in this work that it is the 
nonvanishing P|, i.e., the momentum or the energy associated with the anisotropy change, 
that essentially renders the occurrence of both classical and quantum wormholes possible. 
After all, wormholes seem to be one of the rare means which are under our control (to some 
extent) to challenge the quantum gravity. 
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